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Abstract. Precipitation shifts can have large impacts on hu- 
man society and ecosystems. Many aspects of how inho- 
mogeneous radiative forcings influence precipitation remain 
unclear, however. Here we investigate regional precipitation 
responses to various forcings imposed in different latitude 
bands in a climate model. We find that several regions show 
strong, significant responses to most forcings, but that the 
magnitude and even the sign depends upon the forcing loca- 
tion and type. Aerosol and ozone forcings typically induce 
larger responses than equivalent carbon dioxide (CCL) forc- 
ing, and the influence of remote forcings often outweighs that 
of local forcings. Consistent with this, ozone and especially 
aerosols contribute greatly to precipitation changes over the 
Sahel and South and East Asia in historical simulations, and 
inclusion of aerosols greatly increases the agreement with 
observed trends in these areas, which cannot be attributed to 
either greenhouse gases or natural forcings. Estimates of pre- 
cipitation responses derived from multiplying our Regional 
Precipitation Potentials (RPP; the response per unit forcing 
relationships) by historical forcings typically capture the ac- 
tual response in full transient climate simulations fairly well, 
suggesting that these relationships may provide useful met- 
rics. The strong sensitivity to aerosol and ozone forcing sug- 
gests that although some air quality improvements may un- 
mask greenhouse gas-induced warming, they have large ben- 
efits for reducing regional disruption of the hydrologic cycle. 


1 Introduction 

Precipitation projections vary substantially between models 
in many regions (Meehl et al., 2007). While most of the di- 
vergence may be due to differing representations of physi- 
cal processes, some of the intermodel differences may arise 
from inhomogeneous forcing by aerosols and ozone, as these 


forcings vary greatly across models. Studies suggest that 
hemispherically asymmetric forcing perturbs tropical rainfall 
(Chung and Seinfeld, 2005; Rotstayn and Lohmann, 2002; 
Takemura et al., 2005; Ming and Ramaswamy, 2009), and 
that aerosols affect Sahel precipitation (Ackerley et al., 2011; 
Biasutti and Giannini, 2006; Rotstayn and Lohmann, 2002; 
Kawase et al., 2011; Chang et al., 2011) and the Asian mon- 
soon (Ramanathan et al., 2005; Meehl et al., 2008; Wang et 
al., 2009; Bollasina et al., 2011; Shindell et al., 2012). Hence 
it is important to understand the link between the location and 
type of radiative forcing and the precipitation response. Ad- 
ditionally, virtually all comparisons of the impact of multiple 
forcing agents in emissions trading schemes or evaluation of 
mitigation scenarios use metrics based on global mean forc- 
ing (Shine et al., 2005; Fuglestvedt et al., 2010). Thus it is 
important to know how global mean forcing relates to pre- 
cipitation response for inhomogeneous forcing agents. We 
therefore examined the precipitation response in simulations 
driven by historical individual forcings and by forcing local- 
ized within particular latitude bands. We hope that this initial 
study can provide insight into possible methods for study- 
ing precipitation response to localized forcing and perhaps 
encourage multi-model intercomparison to characterize the 
robustness of these relationships. 

2 Climate modeling and analysis methodology 

We examined transient GISS simulations for 1880-2003 us- 
ing the Climate Model Intercomparison Project (CMIP3) 
version of the GISS-ER coupled ocean-atmosphere model 
with individual forcings by well-mixed greenhouse gases 
(WMGHG), ozone, tropospheric aerosol direct effects 
(ADE), tropospheric aerosol indirect effects (AIE), and all 
forcings (ALL; also including land use, solar and volcanic 
forcings) (Hansen et al., 2007b). For these simulations, we 
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report area-averaged regional trends using the mean of a five- 
member ensemble, where each simulation differs only in its 
initial conditions. Trends are calculated using linear regres- 
sion on the ensemble-mean precipitation at each point, with 
points then averaged together. Uncertainties are assessed us- 
ing an analogous five-member control ensemble of simula- 
tions to quantify internal variability levels, with uncertain- 
ties for a particular regional average based on the standard 
deviation of trends across all the points within a region. In 
other words, each location is essentially treated as an inde- 
pendent measurement of internal variability that goes into 
the regional average (though as neighboring points may be 
related, we in fact assume n/2 degrees of freedom where n 
is the number of points within a region and n/2 roughly ac- 
counts for the modeled autocorrelation between points). 

We concentrate analyses on May through September, 
when absolute precipitation responses are typically great- 
est over Northern Hemisphere (NH) land areas (see Ap- 
pendix A for October to April results). Model results are 
compared with precipitation observations taken from a grid- 
ded dataset that has been optimized for analysis of temporal 
trends (Hulme et al., 1998). Trends are mean linear regres- 
sions over locations with at least 3 of the 5 months having 
data for a particular year and at least 2/3 of years having data. 
Results are robust to reasonable changes in these limits. Un- 
certainties are based on the variability of trends across points 
as in the analysis of modeled trends (again using 1 -sigma 
values and assuming n/2 degrees of freedom where n is the 
number of points within a region). 

We also performed a set of simulations with the identi- 
cal model configuration imposing forcing in discrete lati- 
tude bands (Shindell and Faluvegi, 2009): the tropics (30° S- 
30° N), NH mid-latitudes (30° N-60° N), the Arctic (60° N- 
90° N) and the Southern Hemisphere extratropics (30° S- 
90° S). These bands have generally faster atmospheric mix- 
ing within them than across their boundaries, and hence re- 
gional emissions of short-lived species and their precursors 
have their greatest impact on radiative forcing within these 
areas. Furthermore, the temperature response to forcing oc- 
curs primarily within ~30 degrees in latitude but extends 
very far in longitude (Shindell et al., 2010). Analyses of re- 
sults localized at particular longitudes would also be useful, 
however, especially for precipitation, which could be more 
sensitive to local changes than is surface temperature. 

Equilibrium climate simulations were performed examin- 
ing the response to individual forcings imposed over either a 
single one of the indicated latitude bands or a combination 
of bands. The forcings used were increases in CO 2 and black 
carbon (BC), and decreases in sulfate (representative of re- 
flective aerosols) and tropospheric ozone. As there has been 
very little forcing from tropospheric aerosols or ozone in the 
Southern extratropics, only forcing from CO 2 was imposed 
in that band. Localized CO 2 forcing is as an idealized case 
for illustrative purposes and to aid mechanistic understand- 
ing only, as of course such a forcing is unrealistic. When 
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Fig. 1. Spatial patterns of imposed inhomogeneous radiative forc- 
ings (W m~~ at the tropopause). Forcing were imposed with these 
spatial structures both within individual latitude bands indicated by 
the horizontal lines at 30 and 60 degrees or in multiple bands. 


we imposed forcing in multiple bands to test the linearity 
of the response in comparison with the sum of responses to 
the individual band forcings, we also followed historical pat- 
terns. Hence we examined aerosol forcing from 30° S-90° N 
while for CO 2 we examined the response to forcing at all 
latitudes. All forcings were imposed by scaling present-day 
distributions, so that the aerosol and ozone forcings have a 
spatial structure within the bands comparable to historical 
changes. The spatial pattern of the applied forcing is illus- 
trated in Fig. 1. Spatial variations clearly reflect both the lo- 
cations of emissions and geographic differences in physical 
factors such as the underlying surface albedo. Scaling was 
performed uniformly within the bands, with values chosen 
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Fig. 2. June-July-August average precipitation change (mm day 1 ) per unit positive forcing within the given band (W m -). Forcings are 
imposed only in the areas listed (the historical simulation includes global forcing). Tick marks indicate 30° S, 30° N and 60° N. Regional 
analysis, including statistical significance, is presented in Figs. 3 and 4, for the areas shown in panel (i), which are: India/Bangladesh (71° E 
to 94° E, 15° N to 28° N), Southwest China/SE Asia (98° E to 1 10° E, 11° N to 30° N), Southwestern US (120° W to 103° W, 32° N to 
37° N), Eastern US (95° W to 77° W, 34° N to 44° N), the Pacihc Northwest (129° W to 1 15° W, 42° N to 60° N), and the Sahel (17° W to 
38° E, 9° N to 19° N). 
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to give forcings large enough in amplitude to provide a clear 
response (scaling factors ranged from 0.2 to 2.0, yielding at 
least 0.2 W m -2 local forcing within the band in all cases). 

Simulations were integrated for 120 yr, with analysis us- 
ing area-weighted means over the last 80 yr following sta- 
bilization of the climate. Initial conditions are from a 900- 
yr control run with this model. Significance levels (1-sigma) 
are calculated relative to the unforced variations seen in 80- 
yr segments of that control run. At the regional level, this 
method yields uncertainties that are nearly always about half 
those obtained based on the spatial variation of trends across 
boxes with the regions. A narrower uncertainty range in the 
long equilibrium simulations, which supply many realiza- 
tions of the response to the imposed change, in comparison 
with the uncertainty range calculated for the single realiza- 
tion of transient changes in the historical runs or observations 
is quite sensible. Aerosol perturbation experiments were run 
both with and without a simple parameterization for the AIE 
that allows aerosols to affect cloud lifetime only (Hansen et 
al., 2007b). The response to aerosol forcings applied in in- 
dividual latitude bands tended not to be significantly differ- 
ent if the AIE were included. The Asian response for several 
cases was larger, but not significantly so (see below for fur- 
ther discussion). Hence responses for localized aerosols are 
shown for simulations without AIE to allow comparison with 
historical ADE simulations. 

To compare the effects of the various forcings on a com- 
mon scale, precipitation anomalies are normalized by the 
average radiative forcing applied within the forcing area. 
We use the adjusted radiative forcing (allowing stratospheric 
temperatures only to respond) at the tropopause in all cases 
except for the AIE runs where fixed-SST forcing including 
tropospheric response is used, since forcing without tropo- 
spheric adjustment is not defined for at least portions of the 
AIE (including the cloud lifetime effect simulated here). The 
normalization allows forcings imposed over areas of differ- 
ent size to be sensibly compared. The standard deviation is 
also normalized by this radiative forcing value to maintain 
the model’s true signal-to-noise ratio. 


3 Results 

We first present analysis of the response to individual local- 
ized forcings. We then analyze the response in the historical 
transient simulations, using the insights gained from the in- 
dividual, localized forcing runs to inform our interpretation 
of the historical simulations. 

3.1 Large-scale precipitation response to localized 
forcings 

The spatial patterns of precipitation responses to global CO 2 
increases (Fig. 2a) provides a baseline for our analysis. They 
are almost identical to the response to long-term projections 
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Fig. 3. Ensemble mean June-July-August average precipitation 
change (mm day -1 ) per unit positive forcing (Wm -2 ) for the 
GISS-ER projections driven by A1B scenario concentrations and 
emissions (Hansen et al., 2007a). 


under the A1B scenario (Fig. 3), unsurprisingly as that pro- 
jection is dominated by CO 2 forcing at century timescales. 
In both the GISS model and the mean of the CMIP3 models 
(Meehl et al., 2007), boreal summer precipitation increases 
in a band just north of the equator while generally decreas- 
ing in the subtropics. Precipitation also increases at high 
latitudes. Both these features stem from well-understood 
physics and are robust across climate models (Meehl et al., 
2007; Held and Soden, 2006). The GISS model also repro- 
duces many regional trends that are consistent in sign in > 
90 % of CMIP3 models, including decreased precipitation in 
the Mediterranean, the Pacific Northwest, the Southwestern 
US/Northwest Mexico, and along the North Atlantic storm 
track, suggesting that the precipitation response to globally 
uniform forcing is reasonable. In South and East Asia, the 
Sahel, and the Eastern US, however, the CMIP3 models were 
less consistent in the sign of the precipitation changes. 

The precipitation response pattern for idealized tropical 
CO 2 increases is similar to the global case, though mag- 
nitudes are often smaller (Fig. 2a, d). In contrast, NH 
mid-latitude CO 2 increases induce an opposite response of 
the South Asian monsoon and over Northwestern Mexico 
(Fig. 2a, c) and have little effect on the Southern extratropics. 

Tropical ozone forcing causes a generally similar precip- 
itation response pattern to tropical (or global) CO 2 , though 
there are differences in the tropical Atlantic and in the South- 
ern extratropics. Responses typically have a larger magni- 
tude, especially within the tropics (Fig. 2d, f). In contrast, 
tropical sulfate or BC (Fig. 2h, j) induce quite different re- 
sponses than global or tropical CO 2 in most of the tropics 
(though the extratropical pattern is typically similar). In par- 
ticular, the South Asian monsoon shifts west instead of east, 
and the Sahel dries in response to increased tropical BC. Re- 
sponses are also stronger for equivalent tropical sulfate or BC 
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Fig. 4. Monthly mean precipitation responses (mm day -1 ) to regional and global forcings per unit forcing within the given band (W m -2 ). 
Forcing agent is indicated on the left while the location of forcing is indicated on the right. Individual sets of bars show results for May, 
June, July, August and September (top to bottom). Colored values are statistically significant as compared with the 1-sigma variability in the 
unforced control run (gray areas). 


forcing than for tropical CO 2 or ozone. NH mid-latitude sul- 
fate or BC forcing, however, produces a response similar to 
that produced by NH mid-latitude CO 2 , though again with 
enhanced magnitude (Fig. 2c, g, i). Interestingly, NH mid- 
latitudes are generally more sensitive to tropical than to local 
sulfate or BC, and parts of the tropics are more sensitive to 
NH mid-latitude than to local ozone. The tropics generally 
show high sensitivity to forcing located either in the tropics 


or at NH mid-latitudes, though for CO 2 and aerosols sensi- 
tivity is greater to local forcing. We note that the historical 
response (Fig. 2b) shows large shifts in the tropics, e.g. in 
the location of the Inter Tropical Convergence Zone (ITCZ), 
that seem most consistent with the responses to positive CO 2 
and/or negative sulfate forcing (comparison with historical 
trends is discussed in more detail in Sect. 3.3). 
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Fig. 5. Precipitation sensitivities (left of solid vertical line) and 1900-1998 historical trends (right of solid vertical line) by region. Values are 
May-September averages, given as changes per unit tropopause forcing within the band indicated on the x axis for sensitivities (mm day - 1 
per W m -2 ) and total trend (mm day - 1 ) for historical simulations. Background shading indicates forcing applied in area where precipitation 
response is analyzed. Estimated responses (“Est”; solid symbols) are the sum of the model’s historical forcing times the response per unit 
forcing for each agent and region (using the response to sulfate to represent all reflective aerosols). ADE is tropospheric aerosol direct effects 
and AIE is tropospheric aerosol indirect effects. Obs = observations from (Hulrne et al., 1998). 


3.2 Regional precipitation response to individual 
forcings 

We concentrate on six regions with substantial, coherent bo- 
real summer precipitation anomalies in several experiments 
(see Fig. 2i). Positive forcing at Northern middle or high lat- 
itudes by any agent produces a shift of rainfall away from 
Southwest China/Southeast Asia (SE Asia hereafter) and into 
India/Bangladesh (Figs. 4 and 5). Local tropical BC or sul- 
fate produces an opposite response in India/Bangladesh com- 
pared with tropical CCL or ozone, however. Thus local and 
remote BC and sulfate positive forcing both increase In- 
dia/Bangladesh precipitation, while local and remote CO 2 
and ozone offset one another, making summer precipitation, 
which occurs primarily via the monsoon, highly sensitive to 
aerosol forcing in this region in our model. Note that the re- 
sponses are usually consistent across the months of the wet 


season (Fig. 4). The influence of surface forcing, which can 
affect Asian precipitation (Wang et al., 2009; Ramanathan 
and Carmichael, 2008; Meehl et al., 2008), could account for 
the different response to aerosols versus greenhouse gases. 

SE Asia precipitation responds most strongly to local CO 2 
and ozone, and remote BC. Responses per unit tropopause 
forcing are similar for remote CO 2 , ozone and sulfate and 
typically opposite to responses to local forcing. The response 
to Southern Hemisphere extratropical CO 2 forcing is com- 
paratively small. Per unit surface forcing, responses vary 
greatly across forcing agents in all regions where forcing is 
imposed (Table 1 ; for these and other response areas), sug- 
gesting that change per unit surface forcing is a less useful 
metric for evaluating precipitation response. Hence hereafter 
responses are solely analyzed relative to tropopause forcing. 

For North America, local positive forcing by any agent 
leads to a drier Pacific Northwest (NW) and wetter Eastern 
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Table 1. May-September average precipitation changes per unit forcing (tropopause or surface) within the indicated band (mm day 1 per 
W m -2 ). 



India/ 

India/ 

SW China/ 

SW China/ 

Forcing agent 

Bangladesh, 

Bangladesh, 

SE Asia, 

SE Asia, 

and location 

tropopause 

surface 

tropopause 

surface 


forcing 

forcing 

forcing 

forcing 

60° N-90° N 

BC 

0.15 ±0.04 

-0.22 ±0.05 

-0.09 ±0.03 

0.14 ±0.05 

C0 2 

0.02 ±0.01 

0.04 ±0.03 

-0.02 ±0.01 

-0.04 ±0.03 

Ozone 

0.16 ±0.06 

1.48 ±0.54 

-0.06 ±0.05 

-0.50 ±0.46 

Sulfate 

0.10±0.13 

0. 14 ± 0.18 

-0.05 ±0.11 

-0.06 ±0.15 

30° N-60° N 

BC 

0.41 ±0.04 

-0.21 ±0.02 

-0.20 ±0.03 

0.10 ±0.02 

C0 2 

0.26 ±0.01 

0.67 ± 0.03 

-0.06 ±0.01 

-0.17 ±0.02 

Ozone 

0.18 ± 0.10 

2.92 ± 1.61 

0.03 ±0.09 

0.41 ±1.38 

Sulfate 

0.31 ±0.06 

0.32 ±0.06 

-0.09 ±0.05 

-0.10 ±0.05 

30° S-30° N 

BC 

0.20 ±0.11 

-0.09 ±0.05 

-0.01 ±0.10 

0.01 ±0.04 

C0 2 

-0.16 ±0.03 

-0.78 ±0.13 

0.18 ±0.02 

0.89 ±0.11 

Ozone 

-0.28 ±0.07 

3.42 ±0.90 

0.21 ±0.06 

-2.60 ±0.77 

Sulfate 

0.27 ±0.14 

0.27 ±0.14 

0.14±0.12 

0.14 ± 0.12 

30° S-90° S 

C0 2 

-0.06 ±0.01 

-0.16 ±0.04 

0.03 ±0.01 

0.08 ± 0.03 

30° S-90° N 

BC 

0.78 ±0.23 

-0.44 ±0.13 

-0.09 ±0.19 

0.05 ±0.11 

Sulfate 

0.64 ±0.10 

0.67 ±0.11 

-0.11 ±0.09 

-0.11 ±0.09 

CO 2 all latitudes 

-0.02 ±0.02 

-0.07 ± 0.07 

0.18 ±0.02 

0.66 ±0.06 


US (Figs. 4 and 5). Most tropical forcings lead to similar re- 
sponses, though significance is generally low except for BC, 
and responses sometimes vary across the months of the wet 
season. When forcings are imposed broadly (30° S-90° N or 
globally), BC has a much larger Pacific NW impact than sul- 
fate or CO 2 (per unit forcing), while BC and sulfate have 
greater impacts on the Eastern US than CO 2 . Southwest US 
rainfall is only weakly sensitive to local forcing, but responds 
strongly to tropical aerosol forcing. Responses to tropical and 
broad forcings are the same sign. The strong tropical influ- 
ence on Southwest US precipitation is consistent with pale- 
oclimate studies and sensitivity to tropical ENSO variability 
(Schubert et al., 2004; Meehl and Hu, 2006). 

Sahel rainfall is as or more sensitive to mid-latitude forc- 
ing than to local tropical forcing. Tropical and NH mid- 
latitude BC forcings have opposite impacts (with tropical BC 
showing impacts that vary between the early and late parts 
of the season), and hence the response to broadly imposed 
forcing is greatest for sulfate (and ozone). Substantial Sahel 
drying for the negative forcing caused by increased sulfate 
is consistent with other studies linking observed 1970s and 


1980s drought with NH sulfate (Kawase et al., 2010; Ack- 
erley et al., 2011; Biasutti and Giannini, 2006; Rotstayn and 
Lohmann, 2002; Chang et al., 201 1). 

Since responses are analyzed in terms of precipitation 
change per unit forcing within individual bands, comparison 
of the response to broadly imposed forcings with the sum of 
the responses over the component bands tests the linearity 
of the precipitation changes. With uniform scaling of current 
aerosol distributions, our broadly imposed l.OWm -2 forc- 
ings had roughly 0.6 W m~ 2 tropical forcing and 1.7 W m -2 
Northern extratropical forcing. For the regions examined 
here, responses to broadly imposed CO 2 , BC or sulfate are 
nearly always very near the sum of the responses to the iden- 
tical forcings in the individual bands (Fig. 5; taking into ac- 
count the applied aerosol forcing distributions). The only 
exception is the Eastern US response to sulfate, which is 
slightly larger for the broadly imposed case. 

We call the forcing/response relationships derived from 
these runs Regional Precipitation Potentials (RPPs). These 
provide a way to estimate the regional precipitation re- 
sponse to a given pattern of forcing, analogous to the 
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Regional Temperature Potentials described previously (Shin- 
dell and Faluvegi, 2009, 2010; Voulgarakis and Shindell, 
2010). While there was some dependence on the type of forc- 
ing imposed for the temperature case, the precipitation po- 
tentials vary much more strongly with the forcing agent. We 
examine the applicability of these metrics in the next section. 

3.3 Comparison with historical trends 

We analyzed the response to individual and combined his- 
torical forcings in the transient GISS CMIP3 simulations 
and compared against observations and against estimates ob- 
tained using the response per unit latitude band forcing de- 
scribed above (Fig. 5). The left portion of each panel in the 
figure shows the sensitivity to local forcings as well as three 
cases where forcings were combined across bands (30° S- 
90° N for BC and sulfate, all latitudes for CO 2 ). The right 
portion of each panel shows the response in the historical 
simulations driven by individual and all forcings (open sym- 
bols), the observed regional response (far right symbol), and 
the sum of the RPP (responses per unit local forcing shown 
on the left) times the actual regional forcings imposed in the 
historical simulations (solid symbols). The solid and open 
symbols of the same shape thus show how similar the ac- 
tual response to temporally and spatially varying forcings is 
to the response estimated based on the forcing patterns us- 
ing the RPPs. Comparison of the sensitivities on the left with 
the historical simulation trends on the right helps reveal how 
the various forcing agents in various regions influenced the 
overall trends. Analyses of transient responses used linear re- 
gressions, though for impacts of episodic volcanic forcing 
this may be of limited value. Note that historical WMGHG 
forcing was 2.3 W m -2 , so that WMGHG forcing induced 
changes similar to the C 02 -only case in the historical simu- 
lations as they are typically about 2.3 times larger. 

For regions where observations show substantial trends 
over the 20th century, aerosol forcing is required to cap- 
ture those trends (Fig. 5). In particular, for SE Asia, forc- 
ing by aerosols (including AIE) induced a strong dry- 
ing trend that clearly contributes to the “All forcings” 
response agreeing with observations despite WMGHGs 
causing a large precipitation increase. Over the 1950- 
1998 period, the model’s SE Asia “All forcing” trend is 
—0.14 ±0.07 mm day -1 , in agreement with the observed 
trend of — 0.30 ± 0.21 mm day -1 , again with aerosol-induced 
precipitation trends (—0.50 mm day -1 ; primarily due to di- 
rect effects of local sulfate and remote BC) more than off- 
setting the impact of WMGHGs (0.40mm day -1 ) (consis- 
tent with another recent analysis; Bollasina et al., 201 1). For 
the Sahel, though the model does not fully capture observed 
drying (which is also seen over the past 50 yr; Held et al., 
2005), tropical and NH mid-latitude sulfate increases cause 
substantial drying, while WMGHG, ozone and natural forc- 
ings cause little change or precipitation increases. 


Comparison with the sensitivity studies provides insight 
into the GISS model’s historical simulations (Fig. 2b) in 
other regions as well. Pacific Northwest trends are domi- 
nated by NH mid-latitude BC forcing, the Eastern US shows 
comparable magnitude impacts from WMGHG, ozone and 
aerosols (though BC and sulfate offset one another), while 
for the Southwest US tropical sulfate caused increased pre- 
cipitation, offsetting some of the WMGHG (and BC) induced 
decreases (Figs. 4 and 5). Hence reduction of tropical sulfate 
not only unmasks WMGHG warming, but may also unmask 
WMGHG-induced drying of the Southwest US (though re- 
ducing Sahel drought risk). Decreased precipitation over In- 
dia/Bangladesh, which is consistent with observations over 
the past 50 yr (Wang et al., 2009), is driven largely by aerosol 
forcing, which outweighs the influence of WMGHG, ozone 
or natural forcings. Hence over the next several decades, 
when aerosol forcing is likely to be large, regional precipi- 
tation changes in many parts of the world may be strongly 
influenced by the short-lived climate forcers (as in some pro- 
jections; Levy et al., 2008). 

For ozone changes, the historical response in the tran- 
sient simulations is statistically indistinguishable from an es- 
timate based on the RPPs (a linear sum of the responses 
to tropical and NH mid-latitude forcings times the regional 
historic forcings) (Fig. 5; open and solid diamonds). Note 
that high-latitude ozone forcing was not included in the es- 
timate as historical forcing there stems largely from strato- 
spheric depletion. Similarly, responses to all aerosol direct 
forcings in the transient coupled ocean-atmosphere model 
are consistent with the RPP-based estimates (comparing the 
sum of modeled historical BC and reflective aerosol forc- 
ings times the response per unit forcing for 30° S-90° N BC 
and sulfate forcings) (Fig. 5; open and solid upward point- 
ing triangles). The “All” estimated response is the ADE plus 
ozone estimates added to the actual historic response to other 
forcings (AIE, land-use, solar and volcanic; small forcings 
from stratospheric water vapor and stratospheric ozone are 
included in the All simulation but are neglected in the esti- 
mated response). Non-linearity is seen when all forcings are 
imposed simultaneously relative to the sum of the individual 
components (Fig. 5; open and solid squares) for the Pacific 
NW region, but not for other areas. Hence the RPPs appear 
to provide a useful metric for regional average precipitation 
response in this model, and as with the response to broad 
forcings by a single agent versus the sum of forcings within 
particular bands, responses appear to be generally linear for 
combinations of forcing agents. 

The response in Asia to aerosol forcing when AIE was 
included was larger for several cases relative to simulations 
without AIE (India/Bangladesh response to sulfate, SE Asia 
response to BC), but not significantly so (not shown). An 
enhanced response in those areas is consistent in direction 
with the substantial influence of the AIE in the GISS histori- 
cal simulations (Fig. 5) and with results of another modeling 
study (Bollasina et al., 2011). The AIE-induced changes in 
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the simulations with individual forcing agents are not large 
enough to fully account for that response, however, suggest- 
ing that the net effect may be sensitive to either the mix- 
ture of aerosol types or to types not examined in the indi- 
vidual forcing agent runs. That could include organic carbon 
and/or nitrate aerosols. The impact of regional aerosol in- 
direct cloud effects on precipitation clearly requires a great 
deal more study, especially as the GISS AR4 model only in- 
cluded a highly parameterized cloud lifetime AIE (and semi- 
direct effects) but not cloud albedo effects, and this will be 
the topic of work with the new AR5 model which includes a 
more physically realistic treatment of both types of AIE. 

4 Physical mechanisms 

As noted in Sect. 3.1, forcing by any agent at NH mid- 
latitudes induced a fairly similar pattern of precipitation re- 
sponse in the tropics, with a northward shift of the ITCZ 
during boreal summer. Similar behavior has been reported 
in both other comprehensive climate models and in simpler 
models, and is attributed to thermal forcing and the resultant 
impact on poleward energy fluxes (Chiang and Bitz, 2005; 
Kang et al., 2009). Indeed, such modifications of large scale 
energy fluxes were shown for these simulations in our prior 
analysis of the impact of regional forcings on temperature 
(Shindell and Faluvegi, 2009). Forcing by NH mid-latitude 
ozone is skewed more towards the southern edge of the band 
than the other forcing agents, which may account for its 
stronger impact on the tropics. This result highlights how the 
spatial gradients imposed within the bands and near the lati- 
tude band boundaries may have important impacts on the re- 
sults, and thus merit further study (though we note that the 
quasi-linearity seen in the simulations with multiple band 
forcings suggests these edge effects are not too large for 
aerosols, or at least cancel when realistic perturbations are 
imposed). 

A more detailed examination of the physical mechanisms 
requires analysis of both local and large-scale circulation 
processes. A framework for understanding various contri- 
butions to regional circulation changes in the tropics was 
recently presented by Bollasina et al. (2011). In that the- 
oretical framework, which focuses on the Asia/Pacific sec- 
tor, thermodynamically-induced changes largely follow the 
tropical mean surface temperature anomaly, in accord with 
basic physical understanding (e.g. Held and Soden, 2006). 
The primary thermodynamic impact of warmer tropical mean 
surface temperatures is to slow the tropical zonal circula- 
tion. Dynamic changes (i.e. changes not caused by surface 
temperature or atmospheric moisture anomalies) are induced 
much more effectively by forcing gradients as the circula- 
tion responds to the resulting regional energy imbalances. 
Thus aerosols in the Bollasina et al. (2011) analysis, which 
are primarily cooling agents, by being located primarily in 
the NH act to reduce meridional energy gradients and hence 


substantially slow the merdional circulation while having 
only a modest effect on the zonal flow. The combination of 
aerosol and greenhouse gas forcing reflected the large impact 
of aerosol on meridional circulation and of greenhouse gases 
on zonal circulation in their model. 

We have examined the surface wind and sea-level pres- 
sure (SLP) changes induced by the regional forcings in our 
simulations (Fig. 6). Note that to facilitate comparison of 
the response to various forcing agents, responses to indi- 
vidual forcing agents are all shown per unit positive forc- 
ing (even when the imposed forcing was negative). We find 
that positive tropical BC forcing induces a very strong en- 
hancement of the meridional surface winds over the western 
Indian Ocean, though this is in part reflective of an increas- 
ing anticyclonic flow around the Indian subcontinent. There 
is also a modest enhancement of the equatorial zonal circu- 
lation. Positive tropical sulfate forcing induces a similar but 
weaker meridional response than BC, but causes strong east- 
erly anomalies in the western Pacific and westerly anoma- 
lies in the Atlantic and eastern Pacific sectors. These results 
are somewhat consistent with the framework of Bollasina et 
al. (201 1), especially the zonal flow changes, bearing in mind 
that they examined the impact of an overall negative aerosol 
forcing, but our changes in meridional circulation are in gen- 
eral more closely tied to anomalous flow around localized 
SLP changes. 

Not only do local forcings induce these type of tropical 
responses, but all the NH mid-latitude forcings lead to in- 
creased anticyclonic flow around the Indian subcontinent as 
well. Hence this region is sensitive to both local and remote 
forcing. The response to remote forcing is similar regardless 
of the forcing agent, consistent with thermal forcing, while 
as noted, the response to local forcing is more sensitive to 
BC and sulfate than to carbon dioxide, consistent with dy- 
namic forcing. The enhanced anticyclonic flow increases the 
supply of moisture-laden air flowing from the Indian Ocean 
to the continent during the monsoon, consistent with the pre- 
cipitation responses (Fig. 5). In addition to the Indian Ocean 
sector, there is also a large enhancement of the meridional 
flow off the west coast of Mexico in response to tropical BC. 
This appears to strengthen ascent off the coast, reducing the 
influx of moist marine air into the SW US and accounting for 
the high sensitivity of that region to tropical BC in particu- 
lar (Fig. 5). A somewhat similar response is seen for tropical 
CO 2 forcing, which also strengthens the anticyclonic flow 
here, though the impact of BC is stronger. In contrast, mid- 
latitude forcing has little impact on circulation in the North 
American monsoon area. The enhanced response to BC seen 
in the North American monsoon region is like that seen for 
the South Asian region, though the latter is even more pro- 
nounced, and may reflect the extra energy absorbed within 
the atmosphere by BC. 

Like the tropical responses, extratropical responses are 
also tied to localized circulation changes. The NH mid- 
latitude response to tropical forcing shows anomalous 
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Fig. 6. Sea-level pressure (colors, hPa) and surface wind (vectors) responses per unit positive forcing (top three rows) and in historical simu- 
lations (bottom row; greenhouse forcing is positive, net aerosol forcing negative). Wind changes are normalized to show relative geographic 
variations within a given panel (absolute wind strength changes are closely linked to the magnitude of the SLP anomalies). 


circulation around enhanced North Pacific and North At- 
lantic highs in response to positive forcing by aerosols or 
CO2 and a dipole with positive pressure anomalies over East- 
ern Europe and negative anomalies over the Middle East. NH 
mid-latitude responses over the oceans are similar for posi- 


tive local NH mid-latitude forcing by various agents, and of- 
ten for positive tropical forcing as well, suggesting that in 
these areas positive forcing enhances preexisting circulation 
patterns regardless of the type of forcing or where the forc- 
ing itself is located (although the European sector response 
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is dependent upon forcing location). Responses are typically 
stronger for aerosol forcing than for equivalent CO 2 forc- 
ing, however, leading to the enhanced sensitivity to broadly 
distributed aerosol forcing for the Eastern US, for example 
(Fig. 5). The strong enhancement of the North Pacific high in 
response to BC, which favors southerly rather than easterly 
flow off the coast, is consistent with the reduced precipitation 
reaching the Pacific NW in response to BC forcing (Fig. 5). 

While the tropical responses to BC are enhanced relative 
to other forcing agents, this is typically not the case over NH 
midlatitudes (other than Europe), where the response to BC 
very much resembles the response to sulfate, and hence the 
net aerosol historical impact resembles the sulfate impact (as 
sulfate forcing is stronger than BC’s). Similarly, the historical 
response over the tropical Atlantic does not follow the pattern 
of BC-induced changes but instead matches the response to 
sulfate, especially NH mid-latitude sulfate. 

This suggests that for responses governed primarily by 
thermodynamic forcing, for which all agents cause similar 
responses, the larger radiative forcing by sulfate causes its 
impact to outweigh the impact of BC. This holds for the 
mid-latitudes, and also for western Africa and the Sahel re- 
gion, which is strongly affected by the enhanced equatorial 
north-easterlies induced by negative sulfate forcing (Fig. 6) 
which lead to substantial drying (Fig. 5; consistent with 
thermal forcing). In dynamically-driven monsoon systems, 
BC appears to have an enhanced effectiveness, and it’s im- 
pacts can be as large as those of even a larger sulfate forc- 
ing. Quasi-uniform WMGHG forcing enhances the same NH 
mid-latitude pre-existing circulation patterns as do aerosols, 
but has much less impact on the tropical Atlantic or Indian 
Ocean regions. Ozone forcing, which is much more uniform 
than aerosol forcing (Fig. 1), has broadly similar effects to 
WMGHGs in South Asia, suggesting that it does not trigger 
similar dynamic responses to aerosols. However, its hemi- 
spheric asymmetry still leads to large-scale thermal forcing, 
and it causes a substantial shift in the tropical Atlantic similar 
to that due to aerosols (Fig. 5). 

Thus the meridional circulation anomalies in the Atlantic 
seem to be most closely tied to large-scale (e.g. interhemi- 
spheric) gradients, those in the Indian sector and the North 
American monsoon region are related to both large-scale and 
local forcing gradients, while the anomalies in most other 
regions appear to follow the overall mean forcing magni- 
tude more closely. This suggests that further work explor- 
ing the impact of locating forcing gradients at different lo- 
cations than the bands chosen here would be valuable, es- 
pecially for better understanding Atlantic and Indian sector 
responses (though as noted previously, these sum across re- 
gions fairly linearly, and meridional gradients are compara- 
bly large even for realistic forcings imposed across multiple 
bands for sulfate and BC (Fig. 1)). Intriguingly, the localized 
response over the Arabian Sea is stronger for BC, which has 
very large gradients at these longitudes, while the response 
in the western Pacific is stronger for sulfate, which has its 


own greatest gradients at those longitudes (Fig. 1). Hence for 
dynamical responses, there is a suggestion that the regional 
forcing pattern may be important. 

Our results thus suggest that the gradual shift in aerosol 
and ozone precursor emissions from NH mid-latitudes to 
tropical latitudes that is already underway may have a large 
impact on future precipitation changes in western Africa and 
South Asia, as well as on the Southwest US, while NH mid- 
latitude trends will largely follow the total forcing and be less 
dependent upon its location. Further analyses, including ex- 
amination of water vapor (Seager et al., 2010) and energy 
(Muller and O’Gorman, 2011) budgets as well as the role 
of induced sea-surface temperature changes (Hoerling et al., 
2006) may help further elucidate the underlying mechanisms 
governing the precipitation responses described here. 

5 Discussion and conclusions 

Our regional analyses reveal that precipitation responses to 
aerosols, and, to a lesser extent, to ozone, are often substan- 
tially greater than responses to equivalent CO 2 forcing. This 
result is consistent with initial evidence from the next gener- 
ation of models (CMIP 5th phase) that have shown a much 
stronger sensitivity of South Asian and Sahel precipitation to 
forcing by BC than to equivalent WMGHG forcing (Shin- 
dell et al., 2012) and a strong influence of aerosols on At- 
lantic sector climate trends during the 20th century (Booth 
et al., 2012). Our finding that the heterogeneity of aerosol 
and ozone forcing induces substantial changes in circula- 
tion is also consistent with other analyses (e.g. Rotstayn and 
Fohmann, 2002; Ramanathan et al., 2005; Allen et al., 2012). 
Since the analyses presented here demonstrate that regional 
response patterns depend on the location of forcing, metrics 
based on global mean values do not provide useful informa- 
tion on precipitation impacts for inhomogeneous forcings. 

Our finding that precipitation responses are often linear 
with respect to forcings in different locations and by differ- 
ent agents suggests the possibility of providing useful esti- 
mates of precipitation responses to regional forcing changes, 
however. We show how RPPs that vary with forcing agent 
and location can provide useful information, and in particu- 
lar that the RPP-based estimates of precipitation changes due 
to aerosols and ozone gave values that were generally quite 
consistent with the actual response in full transient climate 
simulations (Fig. 5). This suggests that it might be possible 
to rapidly and easily estimate the precipitation changes in- 
duced by a particular forcing pattern (i.e. without running a 
full climate simulation). Such estimates could provide useful 
guidance for determining which emission scenarios would be 
most useful to explore with comprehensive climate model- 
ing, for example. Future work with higher resolution models 
and multi-model intercomparisons are needed to more fully 
characterize the relationship between localized forcings and 
precipitation response and the variation of the RPP metric 


www.atmos-chem-phys.net/12/6969/2012/ 


Atmos. Chem. Phys., 12, 6969-6982, 2012 



6980 


D. T. Shindell et al.: Precipitation response to regional radiative forcing 


across models, as well as the mechanisms responsible. Fur- 
ther work should also investigate the response to forcings im- 
posed at various longitudes. 

Disagreement among CMIP3 models’ precipitation pro- 
jections is large over North America, the Sahel, and much of 
South and East Asia, where inhomogeneous forcing agents 
have a large impact, but small over Southern Africa, the 
Mediterranean, the North Atlantic, and high latitudes, where 
their impact does not appear to be as enhanced. Differing rep- 
resentations of inhomogeneous forcing may thus contribute 
substantially to divergence in projected precipitation in the 
former regions, consistent with a recent analysis of global 
scale precipitation projections (Pendergrass and Hartmann, 
2012 ). 

An improved understanding of aerosol and ozone influ- 
ences on regional precipitation leads to a more nuanced view 
of the interplay between air quality policies and climate. For 
reflective aerosols in particular, reduced emissions improve 
air quality but augment global warming. This creates a per- 
ception of sulfur controls, for example, having air quality 
benefits but climate disbenefits (Ramanathan and Xu, 2010; 
Raes and Seinfeld, 2009). The response of Sahel rainfall to 
decreased sulfate is a clear example of reduced NH sulfur 
emissions alleviating drought (after the NH emissions first 
caused much of the drought). Thus, although those sulfur 
reductions unmasked greenhouse warming, and hence were 
harmful from a global mean temperature perspective, they 
were extremely beneficial for Sahel precipitation. 

Our results thus imply that policies to address air quality, 
or policies to address climate change that have a substan- 
tial impact on short-lived climate forcers (e.g. fuel switching 
from coal to gas), may have large impacts on regional wa- 
ter supplies that are not proportional to their effect on global 
mean surface temperature. From an air quality perspective, 
any reduction in emissions of aerosol or ozone precursors is 
beneficial. Our analysis suggests that for emissions leading 
to positive forcing, and even for cooling aerosols such as sul- 
fate or organic carbon (OC), reductions may be beneficial for 
climate when considering precipitation in addition to temper- 
ature, since disruption of traditional hydrologic patterns for 
which agriculture and water systems are optimized tends to 
be harmful in the net. This suggests a higher priority should 
be placed on strategies that reduce sulfur and OC than might 
otherwise have been done, as well as reinforcing the evidence 
that reducing BC and ozone precursors is beneficial for mit- 
igating climate change impacts. Our results also imply that 
international cooperation is essential to mitigate hydrologic 
cycle disruptions, as responses are often greater for remote 
forcing than for local forcing. 


Appendix A 

Results for October to April 

We concentrated our analyses on May-September. From 
October-April, modest responses were found in several 
cases in the regions we focus on. The only statistically 
significant response for India/Bangladesh is a slight in- 
crease in response to 30° S-90° N sulfate positive forcing 
(0.04 ±0.02 mm day -1 per W m -2 ). Hence the aerosols ap- 
pear to have a much weaker impact during this season in 
this region than in May-September. For SE Asia, both BC 
and sulfate induce precipitation increases, with tropical forc- 
ings having the largest impacts (tropical ozone has a substan- 
tial impact as well). When imposed broadly (30° S-90° N), 
the increases in responses to BC and sulfate are 0.29 ±0.12 
and 0.10 ± 0.06 mm day -1 per W m -2 , respectively, far more 
than the SE Asia 0.02 ± 0.01 mm day -1 per W m -2 response 
to global COt. In contrast, during May-September, the im- 
pacts of broadly imposed aerosols and CO 2 were of com- 
parable magnitude but opposite sign. The distinctly differ- 
ent behavior for India/Bangladesh and SE Asia during these 
months relative to May-September suggests a response gov- 
erned by different mechanisms than the boreal summer mon- 
soon shift. 

For the Southwest US, October-April responses to 
tropical ozone and sulfate are significant, with both 
causing decreased precipitation (—0.08 ± 0.05 and 
— 0.14±0.10mmday~ 1 per W m -2 , respectively). The 
Eastern US shows increased precipitation in response to 
tropical ozone (0.07 ± 0.03 mm day -1 per W m -2 ) or 30° S- 
90° NBC (0.16 ±0.09 mm day -1 per Wm -2 ), similar to the 
response to global CO 2 (0.10 ± 0.01 mm day -1 per W m -2 ). 
The Pacific NW shows significantly decreased precipitation 
in response to BC imposed broadly or at NH mid-latitudes. 
This response, along with a modest precipitation decrease 
for negative sulfate forcing, appears to dominate the his- 
torical response, outweighing increases in precipitation 
due to WMGHG and ozone forcings. There are almost 
no statistically significant responses for the Sahel during 
October-April, either for individual agents in particular 
bands or for any of the historical simulations. The modest 
responses in October-April relative to May-September in 
most of these areas justify the focus on the latter season. 


Appendix B 

Time dependence of response 

Fooking at transient behavior in first 40 yr of the individual 
localized forcing runs we find that the initial response is gen- 
erally indistinguishable from the equilibrium response. The 
only clear exception is the South Asian monsoon response to 
tropical BC and sulfate, for which the shift in precipitation 
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away from India and into SE Asia in response to positive 
forcing (Fig. 2j, h) takes roughly 7yr to be realized. Our 
results appear to differ somewhat from a prior study (An- 
drews et al., 2010), which found a distinct separation be- 
tween fast atmospherically driven and slow global temper- 
ature change driven precipitation responses to forcing. How- 
ever, their study examined only the global scale, and hence 
included more ocean than land areas, while we focus on land 
areas. As seen in Fig. 2, changes over the ocean are larger 
than those over land, so would dominate at the global scale 
and many of the regional changes we examine are clearly the 
result of large scale shifts in atmospheric circulation rather 
than changes in overall mean rainfall. 
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